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ABSTRACT 

Detailed models of the explosion of a white dwarf, which include self-consistent cal- 
culations of the light curve and spectra, proved a link between observational quantities 
and the underlying explosion model. These calculations assume spherical geometry and 
are based on parameterized descriptions of the burning front. Recently, first multi- 
dimensional calculations for nuclear burning fronts have been performed. Although a 
fully consistent treatment of the burning fronts is beyond the current state of the art, 
these calculations provided a new and better understanding of the physics. Several 
new descriptions for the flame propagation have been proposed by Khokhlov et al. and 
Niemeyer et al.. Using various description for the propagation of a nuclear deflagra- 
tion front, we have studied the influence on the results of previous analyses of Type la 
Supernovae, namely, the nucleosynthesis and structure of the expanding envelope. 

Our calculations are based on a set of delayed detonation models with parameters 
that give a good account of the optical and infrared light curves, and of the spectral 
evolution. In this scenario, the burning front propagates first in a deflagration mode and, 
subsequently, turns into a detonation. The explosions and light curves are calculated 
using a one-dimensional Lagrangian radiation-hydro code, including a detailed nuclear 
network. 

We find that the results of the explosion are rather insensitive to details of the de- 
scription of the deflagration front, even if its speed and the time till the transition to 
detonation differ by almost a factor of two. For a given white dwarf (WD) and a fixed 
transition density, the total production of elements changes by less than 10 %, and the 
distribution in the velocity space changes by less than 7 %. Qualitatively, this insensitiv- 
ity of the final outcome of the explosion on the details of the flame propagation during 
the (slow) deflagration phase can be understood as follows: For plausible variations in 
the speed of the turbulent deflagration, the duration of this phase is several times longer 
than the sound crossing time in the initial WD. Therefore, the energy produced during 
the early nuclear burning can be redistributed over the entire WD causing a slow pre- 
expansion. In this intermediate state the WD is still bound but its binding energy is 
reduced by the amount of nuclear energy. The expansion ratio depends mainly on the 
total amount of burning during the deflagration phase. Consequently, the conditions are 
very similar under which nuclear burning takes place during the subsequent detonation 
phase. In our example, the density and temperature at the the burning front changes 
by less than 3 % , and the expansion velocity changes by less than 10 % . The burning 



conditions are very close to previous calculations which used a constant deflagration ve- 
locity. Based on a comparison with observations, those required low deflagration speeds 
(fti 2 to 3 % of the speed of sound). An exception to the similarity are the innermost 
layers of ^ 0.03 to 0.05 Mq. Still, nuclear burning is in nuclear statistical equihbrium 
(NSE) but the rate of electron capture is larger for the new descriptions for the flame 
propagation. Consequently, the production of very neutron rich isotopes is increased. 
In our example and close to the center, Yg is about 0.44 compared to 0.46 in the model 
with constant deflagration speed. This increases the ^^Ca production by more than a 
factor of 100 to 3.E - 6Mq. 

Conclusions from previous analyses of light curves and spectra on the properties 
of the WD and the explosions will not change and, even with the new descriptions, 
the delayed detonation scenario is consistent with the observations. Namely, the central 
density, results with respect to the chemical structure of the progenitor and the transition 
density from deflagration to detonation do not change. The reason for this similarity is 
the fact that the total amount of burning during the long deflagration phase determines 
the restructuring of the WD prior to the detonation. Therefore, we do not expect that 
the precise, microphysical prescription for the speed of a subsonic burning front has a 
significant effect on the outcome. 

However, on the current level of uncertainties for the burning front, the relation be- 
tween properties of the burning front and of the initial white dwarf cannot be obtained 
from a comparison between observation and theoretical predictions by 1-D models. Mul- 
tidimensional calculations are needed to get an inside into the relations between model 
parameters such as central density and properties of the deflagration front, its rela- 
tion to the transition density between deflagration and detonation, and to make use of 
information on asphericity that is provided by polarization measurements. 

These questions are essential to test, estimate and predict some of the evolutionary 
effects of SNe la and their use as cosmological yard stick. 

Subject headings: Supernovae: general, - nucleosynthesis, abundances - turbulence 
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1. Introduction 

The standard scenario for Type la Supcrnovac con- 
sists of massive carbon-oxygen white dwarfs (WDs) 
with a mass close to the Cliandrasekhar mass which 
accrete through Roche-lobe overflow from an evolved 
companion star (Whelan & Iben 1973; Nomoto & 
Sugimoto 1977). In these accretion models, the ex- 
plosion is triggered by compressional heating. From 
the theoretical standpoint, the key questions are how 
the flame ignites and propagates through the white 
dwarf. Several models within this general scenario 
have been proposed in the past including detonations, 
deflagrations and the delayed detonations, which as- 
sume that the flame starts as a deflagration and turns 
into a detonation later on (Khokhlov 1991, Yamaoka 
et al. 1992, Woosley & Weaver 1994). The latter 
scenario and its variation "pulsating delayed detona- 
tion" , seems to be the most promising one, because, 
from the general properties and the individual light 
curves and spectra, it can account for the majority 
of SNe la events (e.g. Hoflich & Khokhlov 1996, and 
references therein). We note that with the discovery 
of the supersoft X-ray sources, potential progenitors 
have been found (e.g. van den Heuvel et al. 1992; 
Rappaport et al. 1994). 

What we observe as a supernova event is not the 
explosion itself but the light emitted from a rapidly 
expanding envelope produced by the stellar explo- 
sion. As the photosphere recedes, deeper layers of the 
ejecta become visible. A detailed analysis of the light 
curves and spectra gives us the opportunity to deter- 
mine the density, velocity and composition structure 
of the ejecta and provide a direct link to observations. 
A successful application of observational constrains 
requires both accurate early light curves (LC) and 
spectral observations and detailed theoretical models 
which are coupled tightly to the hydrodynamical cal- 
culations (e.g. Harkness 1991, Hoflich, Khokhlov & 
Miiller 1991, Bravo et al. 1996). According to pre- 
vious results, normal bright SNe la can be explained 
by delayed detonation and pulsating delayed detona- 
tion models (e.g. SN 94D, Hoflich 1995). During 
the deflagration phase, the deflagration velocity is 3 
% of the sound speed. In general, a transition from 
deflagration to detonation is required at densities of 
about 2 — 2.5 10'' g cm~^ . Central densities of the 
initial WDs are « 2.0 10^5 cm^^ . Despite their suc- 
cess, the hydrodynamical models are limited by the 
parametrized description of the burning front and the 



ad hoc adjustment of the density at which the defla- 
gration turns into a detonation. 

Recently, significant progress has been made to- 
ward a better understanding of the propagation of 
nuclear burning fronts. First multidimensional hydro- 
dynamic calculations of the deflagration fronts have 
been performed (e.g. Khokhlov 1995, Niemeyer & 
Hillebrandt 1995). A basic, qualitative understand- 
ing of the mechanism which leads to a transition 
from a deflagration to a detonation phase may have 
been achieved (Khokhlov, Oran & Wheeler 1997ab, 
Niemeyer & Woosley 1997). Qualitatively, the re- 
sults agree between different hydrodynamical numer- 
ical simulations but a full description of the deflagra- 
tion in the entire white dwarf and the consistent cal- 
culations of the transition requires high resolution in 
3-D that are beyond the current state of the art. Note 
that the cited references rely on their subgrid mod- 
els and did not resolve the turbulent deflagrations on 
small scales and for appropriate Reynolds numbers. 
Moreover, the transition from a deflagration to a det- 
onation is still not well understood. 

Here, the question is addressed how our results of 
the explosions vary if we use descriptions for the de- 
flagration front which use functional relations derived 
from 3-D calculations. 

1.1. Hydro dynamics 

The explosions are calculated using a one-dimensional 
radiation-hydro code, including nuclear networks (Hoflich 
& Khokhlov 1996, and references therein). This 
code solves the hydrodynamical equations explicitly 
by the piecewisc parabolic method (CoUela & Wood- 
ward 1984) on « 1000 radial zones. The zones near 
the burning front are subdivided to properly track 
its propagation. The code includes the solution of 
the frequency-averaged radiation transport implicitly 
via moment equations, expansion opacities, and a de- 
tailed equation of state. The frequency-averaged vari- 
able Eddington factors and mean opacities are calcu- 
lated by solving the frequency-dependent transport 
equations. About one thousand frequencies (in one 
hundred frequency groups) and about five hundred 
depth points are used. Nuclear burning is taken into 
account using a network of 218 nuclei (see Thiele- 
mann, Nomoto & Hashimoto 1996, Hoflich, Wheeler 
& Thielemann 1998, and references therein). 
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1.2. Description of the Burning Front 

We have considered three cases: 

Case 1) Vhurn = const. Vsound- In our previous in- 
vestigations, const=0.03 has been found to give the 
best fits to observations. This corresponds to the frac- 
tal dimension D= 2 in the description of Woosley & 
Weaver (1994) which suggested I? = 2 to 2.5. 

Case 2 & 3) Here we assumed that Vbum = max{vt,vi) 
where vi and vt are the laminar and turbulent veloci- 
ties, respectively, vi is calculated according Khokhlov 
et al. (1997a). 

Intrinsically, turbulent combustion is a three- di- 
mensional problem. It is driven on large scales by 
the buoyancy of the burning products. The turbu- 
lent cascade penetrates down to very small scales, and 
makes the rate of deflagration independent of the mi- 
crophysics. Turbulent combustion in a uniform grav- 
itational field and static conditions singles out the 
propagation of the flame against gravity g. Chemical 
combustion experiments have been performed in con- 
fined environments, so called Cumbustion chambers. 
These experiments can be reproduced by numerical 
simulations. The propagation speed can be described 
by 
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where ot is the Atwood number, Lf is the character- 
istic length scale, and p+ and p~ are the densities in 
front and behind the front, respectively. However, de- 
spite the success in terrestrial experiments, the basic 
assumptions of both a uniform gravitational field and 
static conditions are violated in the rapidly expanding 
envelopes of SNe la. The main effect of expansion is 
the freeze out of the turbulence on scales L / where the 
turbulent velocity due to Rayleigh Taylor instabilities 
is comparable to the differential expansion velocities 
on those scales, i.e. 



Vt 



yexp 



= Lf/Te; 



Eq.[2] 



Based on this idea, Khokhlov et al. (1997b) sug- 
gested to use the average turbulent velocity (eq. 1), 
use a for uniform, static conditions, and to use the 
mean expansion time scale determined by one dimen- 
sional simulations Texp » dt/d In Rwd- He found for 
the propagation speed of the turbulent burning front 
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time [s] 

Fig. 1. Laminar and turbulent velocities at the 
burning front for models according eq. 1 with C\= 
0.25, 0.20 and 0.15 (top to bottom, m2z02y24i3...1). 
For comparison, VKh gives the velocity of the burning 
front according to Khokhlov (eq. 2, m2z02y24i4). 



Vt = 0M74*^{g Lf) Eq.[3] 

As third case for the description, we followed the 
recipe of Khokhlov but did some modifications by tak- 
ing a, Lf and Texp directly from the hydro at the lo- 
cation of the burning front. Freeze-out was assumed 
when the radius of a mass element has doubled af- 
ter being burned. Ci in equation [1] has been varied. 
Note that a variation in Ci is equivalent to scaling 
the relative length scale for the freeze out. We varied 
Ci in the range to cover a parameter space which in- 
cludes both the descriptions suggested by Khokhlov 
et al. (1997b) and Niemeyer & Woosley (1997). 

2. Results 

The infiuence of the description of the deflagra- 
tion front has been studied at the example of a set 
of delayed detonation model based on the same C/0 
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WD with a mass of 1.39 Mq and a central density 
pc = 2.0 lO^g cm~^. In all cases, a transition density 
ptr of 2.3 10^5 cm~^ has been assumed. The descrip- 
tion of the deflagration front has been varied. The 
deflagration velocity is taken to be 3 % of the speed of 
sound Vsound for m2z02y24i5, and the approximation 
of Khokhlov is used for m2z02y24i4. Eq. (1) has been 
used for models ni2z02y24il-3 with Ci=0.15, 0.20 and 
0.25. M2 indicates that the initial mass of the pro- 
genitor was 2 Mq, with solar metallicity (z02), and 
24 % of Helium. 

In figure 1, the velocity of the burning front is 
shown as a function of time. In general, the speed 
of the burning front is mainly determined by the tur- 
bulent speed but the very early time. As can be ex- 
pected, the transition density is reached later in time 
for smaller Vt because the lower energy production 
per time and, consequently, the slower pre-expansion. 

The final density, velocity and chemical structures 
for the most important elements are given in Figs. 2 
and 3. Overall, the structures are very similar because 
the total energy release depends on the amount of 
the released energy and the initial structure of the 
WD. Even the chemical structure or, more precisely, 
the location of transition between different regimes of 
burning (e.g. from partial to total Si-burning) varies 
between all models by < 7 % in the space of the final 
expansion velocity. The total production of the most 
abundant elements changes by only 4 % and 3 % with 
mean values of 0.60 and 0.16 Mq for M{^^Ni) and 
M{Si), respectively. 




P 110"" g/cm'; day 10] 
V [10' km/sec] 



Fig. 2. — Final density and velocity as a function of 
mass for different models (see text). 

At first, one may expect a rather high sensitivity of 
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Fig. 3. — Same as Fig. 2 but final chemical composi- 
tion as a function of mass. 



the final burning products on details of the descrip- 
tion of the burning front (Hillebrandt, Niemeyer & 
Woosley 1997) but the tests show a low insensitivity 
(Fig. 2 ... 3). This result can be understood by the 
very nature of delayed detonation models. Qualita- 
tively, the low dependence on the details of the flame 
propagation during the (slow) deflagration phase can 
be understood as follows: The production of interme- 
diate mass elements depends on the expansion of the 
outer envelope before the burning front 'arrives'. This 
pre-expansion occurs during the deflagration phase. 
For plausible variations in the speed of the turbu- 
lent deflagration, the duration of this phase is sev- 
eral times longer than the sound crossing time in the 
initial WD. Therefore, the energy produced during 
the early nuclear burning can be redistributed over 
the entire WD causing a smooth lifting/preexpansion. 
In this intermediate state the WD is still boimd but 
its binding energy is reduced by the amount of nu- 
clear energy. The expansion ratio depends mainly on 
the total amount of burning during the deflagration 
phase. Consequently, the conditions are very similar 
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under which nuclear burning takes place during the 
subsequent detonation phase. In Fig. 4, the burning 
conditions just behind the front are given for the two 
extreme models with Ci — 0.15 and Ci = 0.25. The 
durations of the deflagration phase are about 1.7 and 
2.9 seconds, respectively. 



Changes of Burning Condition 




M[MJ 



Fig. 4. — Time at which a mass element is burned, 
mean expansion velocity, density and temperature be- 
hind the burning front for deflagration speeds with 
variable speed according eq. 1 with with various Ci 
(m2z02y24i3/l) and with constant deflagration speed 
(m2z02y24i5). We use arithmetic means over 0.05 

Both models start with laminar deflagrations, i.e. 
the conditions under with nuclear burning takes place 
are similar. After about 3 % of a solar mass has been 
burned, the flame speed is determined by the tur- 
bulent speeds which differ by about 70 %, and the 
densities behind the burning front differ significantly. 
However, in all models, the temperatures are sufR- 
cient high to burn up to NSE and, mainly, ^^Ni is 
produced. After burning of about 0.25 and 0.32 Mq 
for the models with low and high turbulent speeds, 
the transition between the burning as a deflagration 
to a detonation is triggered. Outside « O.4M0, the 
densities and consequently the temperatures are very 
similar. Note that the layers of partial burning are 
located above « 0.6Mq. The expansion rate, rele- 
vant for the time scale for adiabatic cooling, differ 
by 7 % only. This similarity in the burning condi- 



tions explains the insensitivity in the final chemical 
profiles. The pre-expansion is even more similar than 
can be expected from this argument because the hy- 
drodynamical time scales for energy redistribution is 
« 1 second. Thus, the influence of the burning in 
the deflagration mode during the few tenth of a sec- 
ond before turning to a detonation is reduced with 
respect to the pre-expansion of the outer layers. As 
soon as the front turns into a detonation, the remain- 
der of the WD is burned almost 'instantly' or, better, 
on nuclear time scales. 

For comparison, the properties at the burning front 
are given for the description with Vbum/^ sound = 
const. = 0.03 where the constant has been tuned to 
give good agreement between the observed and calcu- 
lated light curves and spectra (e.g Hoflich & Khokhlov 
1996). This description for the burning front is widely 
used in literature and it is consistent with constrains 
from the nucleosynthesis (e.g. Brachwitz et al. 1998, 
Khokhlov 1991, Thielemann et al. 1997, Woosley & 
Weaver 1994 with D=2). 

In the outer layers, the burning conditions are vir- 
tually identical to Ci = 0.15 including the expansion 
rate. Slightly less material is burned during the defla- 
gration phase because the front is faster and the ini- 
tial energy production is larger during the first second 
(Fig. 1.). 

Note that the amount of burning under very high 
densities and, therefore, the production of neutron 
rich isotopes in the central region of the WD depends 
sensitively on the speed of the front which, in the 
more realistic descriptions starts with laminar speed. 
In comparison with previous calculations, this boosts 
the production of very neutron rich isotopes such as 
'^^Ca. In our example and close to the center, Yg 
is about 0.44 compared to 0.46 in the model with 
constant deflagration speed. This increases the '^^Ca 
production from 2 x 10-8Mq to 3 x IO'^'Mq. This in- 
crease is the same for all models which start with the 
laminar deflagration speed. For a systematic study of 
different flame speeds, see the PhD thesis of Brach- 
witz (1999). Note that the production of neutron rich 
isotopes depend on the central density and whether 
the ignition is very close to the center or off-center. 
In principle, this opens a new window for detailed 
analyses of the progenitor and the ignition process in 
SNe la. 
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3. Conclusions 

The final structure of the expanding envelope is 
rather insensitive to the dcrtail(xi (lesc;ription of the 
burning front during tlie deflagration phase. It de- 
pends mainly on global quantities, namely the total 
amount of burning during the deflagration phase and 
the sound speed in the initial white dwarf. Therefore, 
the detailed, spatial structure of the burning front 
cannot be expected to change the result of the explo- 
sion. Findings from previous analyses of light curves 
and spectra on the properties of the WD and the ex- 
plosions will not change and(!) the new description 
of the deflagration front is consistent with the obser- 
vations. The central density, conclusions about the 
chemical structure of the progenitor and the transi- 
tion density from deflagration to detonation do not 
change. Differentially, evolutionary effects and its 
consequences for the observations can be explored by 
parameterized models. 

The validity of the transition density found in pre- 
vious analysis also poses a question on our under- 
standing of the transition from turbulent deflagration 
to detonation which is not well understood. Cur- 
rently, it is attributed to the Zeldovich mechanism 
which is based on the adiabatic mixing of burned and 
unburned material where the entropy increase in the 
unburned fuel triggers the explosion (Khokhlov et al. 
1997b, Niemeyer et al. 1997). The relative change 
over the entire range of parameterizations discussed 
here corresponds to a change in the transition den- 
sity of ~ 10%. This leaves us with the puzzle why 
the theoretical estimates for the transition density are 
lower by a factor of about 0.7 (Khokhlov et al. 1997b) 
and 0.4 (Niemeyer et al. 1997). For reasons, we can 
only speculate. The diflferences between the latter val- 
ues may indicate the size of the uncertainties in our 
understanding of this transition process. Other, not 
yet considered microscopic effects may be involved. 
Another possibility may be that we measure two dif- 
ferent things. To derive the model parameters from 
the observations, we measure the mean density at the 
burning front when the transition occurs whereas the 
theoretical considerations provide information on the 
location where the transition occurs. Maybe the frag- 
mentation of the burning front causes that the transi- 
tion occurs somewhat ahead of the mean front. If this 
interpretation is correct, this may indicate that the 
detonation is started at about 10 to 20 % (in radius) 
ahead of the mean front. In either CcLS6, cl clarification 



needs further investigations and, certainly, will pro- 
vide new inside into the properties of nuclear burning 
fronts. 

We note that the nucleosynthesis, optical and IR 
spectra and light curve seem to require that, after an 
initial phase of slow burning, the front moves with 
velocities close to the speed of sound to keep up with 
the expanding outer layers and to burn Carbon under 
low density conditions. The latter is indicated by a 
strong line at 1.05 yum which can be attributed to Mg 
II and indicates expansion velocities in excess 15,000 
... 16,000 km/sec. In general, unburned Carbon is 
restricted to velocities above 20,000 km/sec (Fisher 
et al. 1997, Hoflich 1995, Hoflich & Khokhlov 1996, 
Khokhlov 1991, Nugent et al. 1997, Thielemann et al. 
1996, Wheeler et al. 1998). Both the high velocities 
in Mg and C are consistent with delayed detonation 
models but inconsistent with the classical deflagration 
model W7 (Nomoto et al. 1984). Extended mixing 
of the inner, Ni-rich layers can be excluded from the 
spectra (Hoflich, Wheeler & Thielemann 1998). Mix- 
ing of the outer layers may produce Si lines at high 
velocity (Nugent et al. 1997) but does not solve the 
problems with Mg and C. Although the evidence fa- 
vors the delayed detonation scenario in the framework 
of 1-D models, this does not necessarily imply that 
a transition from detonation occurs. Alternatively, 
the flame may propagate as a very fast deflagration 
wave (Hillebrandt, private communication). In this 
context, the transition density may be attributed to 
the region where the turbulent flame is not any more 
driven by Rayleigh- Taylor instabilities. 

Finally, we have also to stress the limits of 1-D 
models with a parameterized description of the de- 
flagration front and the need for multidimensional 
hydrodynamical calculations and radiation transport. 
These limits are critical to understand details of the 
explosion and to use SNe la as distance indicators on 
cosmological scales (e.g. Schmidt et al. 1998, Perl- 
mutter et al. 1999). Observations clearly showed a 
strong relation between the maximum brightness and 
the decline relation (Phillips 1993). From theoretical 
models, the amount of radioactive nickel produced 
during the explosion of a massive white dwarf has 
been identified as the basic quantity which provides 
this relation. The relation is independent from details 
of the model, however, the amount of '"'^iVi actually 
produced depends on a combination of free param- 
eters in the models such as central density and the 
chemical composition of the WD, and the propaga- 
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tion of the burning front. If these parameters are 
varied independently within the hmits indicated by 
individual fits to observations (Hoflich & Khokhlov 
1996); we expect a spread around the mean maxi- 
mum brightness/decline relation of « 0.4™ which is 
consistent with the spread based on the CTIO data 
published by Hamuy et al. (1996). Recent redefi- 
nitions of the statistical methods and new observa- 
tions suggest a much tighter relation with a spread of 
w 0.12" (Schmidt et. al. 1998). This narrow spread 
cannot be understood in the context of the parameter 
range allowed by current analyses. This tight relation 
may hint of an underlying coupling of the progeni- 
tor, the accretion rates and the propagation of the 
burning front. However, on the current level of un- 
certainties, the relation between the free parameters 
such as transition density and initial central density 
cannot be deduced from a comparison between ob- 
servation and theoretical predictions by 1-D models. 
Multidimensional calculations are needed to test e.g. 
the relation between chemical composition or central 
density of the white dwarf and the properties of the 
deflagration front and its relation to the transition 
density between deflagration and detonation. 

In general, the chemical signature of the deflagra- 
tion phase is 'wiped out' by the detonation front but 
a small trace may remain at the interfaces between 
complete and incomplete Si burning (Thielemann et 
al. 1996). Although not conclusive, recent polariza- 
tion measurement may already hint toward the exis- 
tence of such inhomogeneities in the chemical struc- 
ture (Wang, Wheeler & Hoflich 1998). 
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